Slow nucleation rates in Chain Inflation with QCD Axions or Monodromy 
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The previous proposal (by two of us) of chain inflation with the QCD axion is shown to faiL The 
proposal involved a series of fast tunneling events, yet here it is shown that tunneling is too slow. 
We calculate the bubble nucleation rates for phase transitions in the thick wall limit, approximating 
the barrier by a triangle. A similar problem arises in realization of chain inflation in the string 
landscape that uses series of minima along the monodromy staircase around the conifold point. The 
basic problem is that the minima of the potential are too far apart to allow rapid enough tunneling 
in these two models. We entertain the possibility of overcoming this problem by modifying the 
gravity sector to a Brans-Dicke theory. However, one would need extremely small values for the 
Brans-Dicke parameter. Many successful alternatives exist, including other "axions" (with mass 
scales not set by QCD) or potentials with comparable heights and widths that do not suffer from 
the problem of slow tunneling and provide successful candidates for chain inflation. 

PACS numbers: 



I. INTRODUCTION 

In the old picture of inflation, originally suggested by 
A. Guth 1], the universe is trapped in a meta-stable 
false vacuum whose energy causes the universe to expand 
quasi-exponentially. As the universe expands, bubbles of 
true vacuum gradually form inside the sea of false vac- 
uum. The transition is completed by the coalescence of 
bubbles, releasing the latent heat of transition and re- 
covering from the supercooling. To have the universe ex- 
pand sufficiently to solve the problems of standard big 
bang cosmology, the phase transition should be quite 
slow. This would prevent the bubbles from percolation 
and thus the universe would never recover from the in- 
flationary stage. 

Chain inflation resolves the problems of old infla- 
tion by assuming that more than one stage of inflation 
is responsible for solving the problems of Big Bang cos- 
mology. Instead, the universe tunnels rapidly through a 
series of ever lower energy vacua. The failure of old infla- 
tion is avoided because the bubbles of true vacuum are 
able to percolate at each step since the phase transition 
is fairly rapid. The nucleation rate per unit four volume 
has the form 0]: 



Ae 



-Se 



(I.l) 



where S'b is the Euclidean action for the bounce solution 
extrapolating between false and true vacua and A is a 
determinant factor which is generally of order the 
quartic power of the energy scale of the phase transition. 
For a flrst order phase transition, with Einstein gravity, it 
has been shown that the probability of a point remaining 
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in the false vacuum is given by 



47r 

p{t) ^ exp(-— /377i), 



where (3 is deflned by 



/9 



(1.2) 



(1.3) 



Writing Eq. (|I.2p as exp(— i/r), the lifetime of the field 
in the false vacuum is given by: 



_ 3 

The number of e-foldings for the tunneling event is 

3 H-^ 



(1.4) 
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Hdt ~ Ht 
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(3 has to be greater than the critical value, /3c, where 



(3c = 9/4^, 



(1.6) 



to achieve percolation and thermalization. This corre- 
sponds to having an upper bound on the number of e- 
foldings that that can be obtained in each stage of infla- 
tion: 



1 

X < Ac = 



(1.7) 



The approximate number of e-foldings which is re- 
quired to solve the problems of standard cosmology de- 
pends on the energy scale of inflation, M, and the re- 
heating temperature: 

iVe = 45 + ln(3000) + ^ ln(Afi4) + ^ In(rio), (1.8) 

where M = MmIO^* GeV, Trh = TiolOi° GeV. For 
a model with energy scale and reheating temperature 
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around the GUT scale, about 60 e-foldings are required, 
whereas for the one in which these parameters take values 
around the QCD scale, minimum number of e-foldings is 
reduced to 25. This corresponds to having at least around 
200 phase transitions for the GUT scale chain inflation 
and around 100 ones for the QCD scale chain inflation. 

A study of density perturbations from chain inflation 
has been performed in [7, 8], which found that the right 
amount of perturbations to match data can be generated. 
Previous studies can be found in [1, [l^, [HI . 

Chain inflation can arise in many contexts. However, 
we wish to make clear at the outset that many other 
versions of chain inflation may be very successful. For 
example, Freese, Liu, and Spolyar tl3| examined chain 
inflation due to four-forms in string theory where there 
can be large numbers of potential minima connected by 
tunneling. We note that the present result does not apply 
to the case of four-form inflation. 

In this paper we focus specifically on two variants of 
chain inflation and show that the model fails for these 
cases: i) chain inflation with the QCD axion fis'l and 
ii) monodromy chain inflation 0, The basic prob- 
lem is that, for these two specific examples, the width of 
the potential is so large that tunneling never takes place. 
Individual minima are (meta)stable for extremely long 
timescales. We use a triangle potential to approximate 
the real potentials in obtaining these results. However, 
if one considers "axions" at any scale, not tied to QCD, 
then chain inflation can work for the scales discussed here 
in the paper. In addition, for potentials with comparable 
heights and widths, the problem can be avoided as well. 
Chain inflation fails for some select models but can easily 
work in others. 

First we will review QCD axion inflation and remind 
the reader how the necessity to avoid the thin wall limit 
arose in the original paper [31 . Then we calculate the 
nucleation rate for the thick wall limit, approximating the 
potential with triangles. Next we will relax the restric- 
tion of taking the inflaton to be the QCD axion. If the 
inflaton is some other axion with different mass scales, 
we can set some constraints on the height and width of 
the bumps that separate two adjacent minima, in order 
to have critical nucleation rates. These constraints could 
be used in the future to identify proper candidates for 
realizing chain inflation. Finally, inspired by extended 
inflation [l6!|, we investigate the possibility of overcom- 
ing the difflculties of QCD axion inflation by modifying 
the gravity sector of the theory. 



II. CHAIN INFLATION WITH THE QCD 
AXION 



provide a tilt to the potential. Here we discuss this model 
in detail and examine its pitfalls as a proposed inflation 
model. 

The QCD axion is the Goldstone boson of the bro- 
ken U{l)pQ Peccei-Quinn symmetry which was postu- 
lated to solve the strong CP problem of strong inter- 
actions llBl- As was shown by Sikivie [l9|, when 
this global symmetry is broken, a discrete Zjv subgroup 
remains unbroken, where N is the number of quark fla- 
vors that rotate under ?7(1)pq. The complete form of 
the axion potential depends on non-perturbative effects. 
However, for definitcness, we will focus on the "invis- 
ible axion" model of Dine-Fischler-Srednicki-Zhitnitsky 
[20l . [2]| , where the axion is identified as the phase of a 
complex SU{2) x U{1) singlet scalar, a = ex.p{ia/v). 

Below the Peccei-Quinn symmetry breaking scale, 
the axion potential takes the following form: 



V{a) = Vo 



1 — cos 



Na 

V 



(II.l) 



Here v ~ faN, where fa is the decay constant of the 
axion; it has to be between 10^ GeV and 10^^ GeV to 
satisfy the constraints from stellar evolution [l^, and 
cosmology [24]. In addition, Vq = m^f^ is on the or- 
der of the QCD scale, which we take to be around 220 
MeV. The axion mass, ma, is equal to 2iV j^/ttTO^/w, 
where z = riiu/md — 0.56 (but see (25j). and m^r = 135 
MeV and /tt = 93 MeV are respectively the mass and 
decay constant of the pion. The energy scale of in- 
flation is around the QCD scale, and thus, as pointed 
out above, the required number of minima has to be 
taken to be larger than 75 to obtain a sufficient num- 
ber of e-foldings. This requires the introduction of ad- 
ditional heavy fermions beyond the usual quark and lep- 
tons. They carry color charge as well as PQ charge and 
thus contribute to the QCD anomaly 19]. To be definite 
hereafter we take A'^ = 100. 

The J7(l)pQ symmetry is softly broken by adding a 
term of the form fjPa + h.c. to the Lagrangian. This in 
turn adds a term of the form rj cos(a/v + 7) to the axion 
potential, where rj and 7 are real parameters. (Note that 
7 misaligns the QCD and soft breaking minima.) Intro- 
ducing a new variable 9 = a/ fa, the combined potential 
for the axion takes the form: 



U(0) = Uo(l-cos0)-r;cos(- 



-7)• 



(H.2) 



There is a limit on 7 from the electric dipole moment 
(EDM) of the neutron, namely, from the fact that the 
minimum of the combined potential should not be shifted 
from zero: 



Two of us [13| previously proposed that the QCD ax- 
ion might be the inflaton if it were supplemented with i) 
a large number of minima, corresponding to new heavy 
fermions participating in the chiral anomaly and ii) soft 
symmetry breaking of the Peccei-Quinn symmetry to 



EDM 



< 6 X 10 
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(II.3) 



We should impose the EDM bounds at the bottom of the 
potential, which is the end point of tunneling. There one 



can show that A9 



Tj Sin 7 
'VoW 



I13II. This constraint can be 



3 



V{(l>) 




FIG. 1: To calculate the nucleation rate we have approxi- 
mated the potential with a triangle. 



satisfied in many ways as was discussed in [13|. Here, for 
simplicity, we set 7 = from the outset. 

In the thin wall limit, the Euclidean action is given by 



Se = 

where Si is given by 
Jo 

Thus, we have: 



2e3 



1 — cos 



la 



da. 



Se = 5x 10^ 



(II.4) 



(II.5) 



(II.6) 



For the parameters of the invisible axion [20], Se ^ 1 
and thus tunneling is suppressed. In it was specu- 
lated that being away from the thin wall limit and having 
e < Voi can cure this problem. However as we will show 
below, the problem of suppressed tunneling persists even 
in the thick wall limit. 



III. THICK WALL LIMIT 

To compute the nucleation rate, we will use the result 
of [26j . where a general one dimensional potential could 
be approximated by a potential with triangular form (see 
Fig. [T|), in which case the equation of motion for the 
bounce can be explicitly solved. Note that ((/)_|_, V+) and 
(0_ , V- ) are the values of the field and potential energies 



at the false and true vacua, respectively, and {(j)T, Vt) is 
the value of the field and potential energy at the maxi- 
mum of the potential. The resulting expression for the 
Euclidean action takes on two forms, depending on the 
condition [35j]: 



AV+ 



1/2 



> 



2A0_ 



(III.1) 



Here we have defined the quantities A(/)± = ±(</)t — 4>±)i 
AV± = {Vt - V±) and A± = Ay±/A0±. If the above 
condition is satisfied, Se takes the form: 



Se = 



327r^ 



1 



A0t 
AV+ 



(111.2) 



where c = A_/A-|_. Otherwise, if condition pil.ip is not 
met, the Euclidean action is instead given by: 



Se = ^^^^+^Rt {-13', 



(in.3) 




Rt 



1 f0l + c0^ 

2 U/3--/34 



(in.4) 



In order to connect the parameters of the axion poten- 
tial to that of the triangle approximation, we first note 
that e, the energy difference between two consecutive 
minima, is proportional to the soft-breaking parameter 
•q. However we cannot make rj arbitrarily large because 
once 77 becomes equal to NVq the barrier between two 
vacua disappears and the field simply rolls down the po- 
tential rather than tunneling. Thus, if N happens to be 
close to the critical value, Nc ~ 100, for which one would 
obtain enough e-foldings, the model will not be able to 
inflate enough if 77 is larger that NcYq. Hereafter, we 
will take 77 = [Nc — l)Vb since it gives the most rapid 
tunneling rate without going over to rolling. 

We note that, in the middle of the axion potential, 
the condition (|III.1|) is satisfied and hence we may use 
Eq. (|III.2[) to compute the instanton action. For fa = 
10^2 GeV (or 10^ GeV), we see that Se is of order lO^^ (or 
10'^'*) in the middle of the potential. Thus the tunneling 
is highly suppressed. At the bottom of the potential, 
where e gets smaller, condition pil.l[) is not met and we 
have to use equation pil.3p to calculate the Euclidean 
action. Here the situation gets worse; for fa — 10^^ GeV 
(or 10^ GeV), the Euclidean action takes the value 10^'^ 
(or 10*^). The minute value of the nucleation rate can be 
understood intuitively by noting that the vacua are very 
far apart and hence their wave functions cannot overlap. 

Somewhere in the limit of a broad barrier, the Coleman 
transition [sf ceases to exist. Then the decay instead goes 
through the Hawking-Moss instanton [23|, which physi- 
cally describes a thermal jump of the field onto the top of 
the barrier. However, to have Hawking-Moss instanton. 
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we must have ^ ^ 2H^ . This condition is not satis- 
fied in the case of the QCD axion, for which ~ ^o//a j 
as fa ^ M-p\. Thus we conclude that the QCD axion, as 
suggested in , is not able to realize chain inflation. 

"Axions" not tied to the QCD scale are also ubiquitous 
in string theory [28l | . If we relax the requirement that we 
are dealing with the QCD axion here, and consider "ax- 
ions" of arbitrary mass scale, successful chain inflation 
models can be found that avoid the problem of insuffi- 
cient tunneling. For a potential of the form (|II.2[) with 
arbitrary barrier height Vq and decay constant /„, the 
instanton action Se can become of 0(1), if the following 
constraint is satisfied: 



fa < lO-'V, 



(III.5) 



It would be worthwhile to investigate the realm of string 
theory axions to see if one can find suitable candidates 
that satisfy (jIII.5[) to realize the idea of chain inflation. 

A similar problem of suppressed nucleation rates arises 
in a particular attempt to realize chain inflation in the 
string landscape [tJ . In [l5| , the generic distance between 
minima in the monodromy staircase is obtained to be 
large, O(Mpi). As requires a shorter inter-minimum 
distance, this will again lead to suppressed nucleation 
rates. 



IV. EXTENDING CHAIN INFLATION 

In the remainder of this paper, we suggest a way of 
avoiding the problem of suppressed tunneling rates of 
QCD and monodromy chain inflation by modifying the 
gravity sector of the theory. The remedy is motivated 
by how extended inflation [l^| amends the graceful exit 
problem of old inflation. 

Several approaches have been taken to solve the grace- 
ful exit problem of old inflation. Adams and Freese [l^ 
suggested a time-dependent nucleation rate, so that the 
initial value of (3 deflned in Eq. p.3p was very small (and 
the universe inflates) but the flnal value of f3 becomes 
very large (so that the bubbles of true vacuum percolate 
and the phase transition completes) ; their model is called 
double field infiation. 

Around the same time. La and Steinhardt [iGj made an 
alternative suggestion, namely extended inflation, where 
the gravity part of the action takes a Brans-Dicke (BD) 
form 30] , 



BD 



(IV.l) 



where lo is known as the Brans-Dicke parameter. In the 
above action, general relativity is recovered in the limit of 
uj —^ oo and present day tests of general relativity require 
uj > 500. Modification of the gravity sector as above will 
change the expansion from de-Sitter to p^ : 



ait) = I 1 + ^ 
' a 



(IV.2) 



Here = [3 + 2uj){5 + 6u;)/12 and = 87rpo/3$(0), 
where $(0) is mpj in the beginning of inflation at t = 
and po — V{Q) is the value of cosmological constant in the 
metastable vacuum. For t < t^s, where ids — Q^/x, the 
scale factor expands exponentially. However after that 
the behavior of the scale factor changes to power-law, 
f^, where q = uj + 1/2. Change in the evolution of the 
background modifies the rate at which bubble nucleation 
converts false vacuum to true vacuum. 

As shown in the probability of a point remaining 
in the false vacuum is given by 



p{t) = exp 



47r 



dt'T{t')a^{t')) — 



3 \Jt' ait"), 

(IV.3) 

For de-Sitter expansion, this reduces to p.2p . However, 
for a power-law background, the probability to remain in 
false vacuum phase becomes: 



pit) = exp [--I3ujy' 



(IV.4) 



where y = x^/^ 0- -^s pit) is decreasing much faster 
than the volume, the universe will soon exit inflation and 
graceful exit is achieved. However, for the case of ex- 
tended inflation |26|, it was understood that the bubbles 
that form in the beginning of inflation would be stretched 
to cosmological scales by inflation. This constraint forced 
Lu to be less than 20. With this value of w, requiring the 
scalar spectral index to satisfy constraints from WMAP 
is impossible. Resolutions that are based on making lo 
time-dependent are also very contrived [3l|. Nonethe- 
less, we may apply this idea to chain inflation to see if 
we can obtain a sensible, even though contrived, solution 
for the small nucleation rate of QCD and monodromy 
chain inflation. 

We first calculate tp, the time it takes for the universe 
to complete the power-law phase. The volume in the false 
vacuum phase is 



(IV.5) 



where pit) is given in Eq. (jIV.4[) . The volume in the false 
vacuum phase starts to contract at time 



tc = 



1/4 



X 



(IV.6) 



At tc, the exponential factor in p(t) starts to takes over 
the behavior of Vf and thereafter it starts to decrease 
exponentially, with an exponent equal to — (t/ic) ■ Only 
within tp ~ few x tc, the volume in the false vacuum 
phase becomes an infinitesimal fraction of the total vol- 
ume. Therefore the coalescence of true vacuum bubbles 
and recovery from the false vacuum phase is achieved 
within few x tc- The number of e- foldings obtained in 
the power-law phase is then: 



Np~qln 



tdS 



In 



47r/3a4 



(IV.7) 
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This should be compared with the number of e-foldings 
obtained in the de-Sitter stage: 

Nds = a. (IV.8) 

Note that Brans-Dicke theories with w ~ 1 arise nat- 
uraUy in compactifications of sting theory and other 
Kaluza-Klein theories [l^l- For uj — 1, we have 

Np ~ Se (IV.9) 
Nds =^ 2.5 (IV. 10) 

where S'^ is the Euchdean action for the bounce solution 
which is obtained in the last section. For QCD inflation, 
at the bottom of the potential, Se was obtained to be of 
order 10^'^ for fa — 10^^ GeV. Therefore the number of e- 
foldings obtained in the power-law phase of the first min- 
imum is much more than enough and can solve the prob- 
lems of standard Big Bang cosmology by itself. There is 
no need for the other minima in chain inflation. One can 
circumvent the long power-law phase by adjusting the 
value of LU to be of order (3 ~ exp(— S's/S) ~ exp(— 10^^). 
For such values of lo the universe expands one e-fold 
in the de-Sitter stage and recovers from the power-law 
phase in a small fraction of an e-fold. Although such 
small values of lo are in disagreement with solar tests of 
general gravity, it may be possible to make the parame- 
ter to field-dependent (similar to hyper-extended inflation 
[1^), so that it takes such small values during inflation 
but is then driven toward experimentally-viable values 
afterwards. 

In summary, we have identified a problem of too slow 
nucleation rates in chain inflation in those models with 



minima separated by great distances in field space. The 
two cases we studied were chain inflation with the QCD 
axion [l^ and monodromy Q chain inflation. To ob- 
tain sufficiently rapid tunneling, one could look for other 
"axionic" candidates that satisfy constraints like (jlll.sp . 
We have also suggested another option of resolving the 
problem by extending the gravity sector of the theory to 
Brans-Dicke theory. This solution is motivated by the 
mechanism whereby extended inflation solves the grace- 
ful exit problem of old inflation (which also arises due to 
the small nucleation rate of the inflaton). Such a reso- 
lution requires an extremely small value for the Brans- 
Dicke parameter, w, which must therefore become time 
dependent. It is interesting to see if string theory can 
come up with such values for lo. Alternatively, many 
models of chain inflation which have smaller distances 
between vacua (or are not described by a potential at 
all) easily have rapid tunneling and do not suffer from 
the problems discussed in this paper at all. 
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